KltS 

'article Physics 



YITP-05-59 
RITS-PP-006 



Graviton emission from a higher-dimensional black hole 

A. S. Cornell, 1 ^ Wade Naylor, 2 'Q and Misao Sasaki 1 '© 

1 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
2 Department of Physics, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan 

(Dated: October 2, 2005) 

We discuss the graviton absorption probability (greybody factor) and the cross-section of a higher- 
dimensional Schwarzschild black hole (BH). We are motivated by the suggestion that a great many 
BHs may be produced at the LHC and bearing this fact in mind, for simplicity, we shall investigate 
the intermediate energy regime for a static Schwarzschild BH. That is, for (2M) 1 ^ n ~ 1 ^ ~ 1, where 
M is the mass of the black hole and u> is the energy of the emitted gravitons in (2 + n)-dimensions. 
To find easily tractable solutions we work in the limit / ^> 1, where / is the angular momentum 
\^ ' quantum number of the graviton. 
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co ' i. introduction 

h ; 

Much discussion has recently focused on the emission rates of TeV BHs as motivated by the proposition that the 
^ ■ quantum gravity scale can be brought down to as low as a TeV in some higher-dimensional BH models [J, 0, Q • Most 
ON ' cases have focused on the low-energy emission of scalar or spinor fields from higher-dimensional Schwarzschild and 
slowly rotating Kerr BHs, for a review see reference Q. Numerical methods also allow us to evaluate the emission in 
the full energy range [^ . 

The importance of bulk emission of gravitons, as well as possible recoil effects, was highlighted recently in references 
0, (also see references 0, H). However, little attention has been paid to the actual emission of gravitons from 
^) \ a BH (though there has been some work relating to the quasi-normal modes (QNMs) of higher-dimensional BHs 
[ToL fllL IT2L PHI ) and furthermore it is useful to have analytic expressions for the cross-sections etc. not just in the 
low-energy or high-energy classical regime. In this article we shall discuss BH cross-sections for gravitons in what 
we shall call the intermediate energy regime for the variable e = 2Muj n ~ 1 , that is where e ~ 1. In this regime the 
energy of the particle is near the peak of the potential barrier in the associated scattering problem. The classical 
cross-section is reproduced in the high-energy limit, £ > 1, as we shall discuss later. 

In this article we shall focus on the static Schwarzschild BH, though in some models a rotating BH does not 
necessarily spin-down to zero, but evolves to a non-zero angular momentum [lJ,LL5|. However, as yet the gravitational 
perturbation equations for a higher-dimensional rotating BH are unknown. Thus, we shall investigate the intermediate 
energy regime for the static case, hoping that it may encode some of the properties of the rotating case. Indeed, 
graviton super-radiance for a large rotation parameter is expected in the intermediate energy regime. To this end only 
recently has the case of spin-zero emission for the high-energy and high angular momentum regime been studied 
llT^ , where numerical methods were employed. Note that scalar emission for low-energies and low angular momentum 
has been studied in references pH . . 

The particular background we shall investigate is that of a static Schwarzschild BH in (2 + n)-dimensions. To 
perform the calculation we shall recall some results recently derived for the gravitational perturbations of a higher- 
dimensional maximally symmetric BH [2(3. As mentioned earlier, although it seems likely that TeV BHs will probably 
be highly rotating when produced 15] there is not yet a method for separating the gravitational perturbations on a 
higher-dimensional Kerr background. Thus, we shall content ourselves with the static Schwarzschild case for now. 

Gravitational perturbations in general separate into scalar (polar) and vector (axial) perturbations. In more than 
four-dimensions there is an extra degree of freedom corresponding to tensor perturbations. As shown in reference [2(3 



o 
o 



* Email: alanc@yukawa.kyoto-u.ac.jp 

t Email: naylor@se.ritsumei.ac.jp 

^ Email: misao@yukawa.kyoto-u.ac.jp 



2 



we have for the scalar perturbations; 1 



dr V dr 

where 



/i(/?'l+V s * = ***, (1.1) 



and 



16H[ra + -|n(n + 1)(1 — /)J 2 

with 

#(r) = (n 4 (n + 1) 2 (1 - ff + n(n + 1) [4(2n 2 - 3n + 4)m + n(n - 2)(n - 4)(n + 1)] (1 - f) 2 

-Yin [(n - 4)m + n(n + l)(n - 2)] m(l - /) + 16m 3 + 4n(n + 2)m 2 ) . (1.4) 

In the above 

rn = l s (ls + n — 1) — n, ^ = 2,3,... (1.5) 

and as discussed in reference [2(J the mode I = corresponds to adding a small mass to the BH; the mode I = 1 is 
a pure g au g e mode. Note that the true mass of the BH, Mbh, as measured on the brane, is the same as that in the 
bulk [23, [24J and hence; 

nA n M 

MBH = 8^0^ ' (L6) 

where 

27r (n+l)/2 

A = r[(n+l)/2] (L7) 

is the area of a unit n-sphere, G n +2 is the (2 + n)-dimensional Newton constant, and c is the speed of light. In what 
follows we shall set G n +2 = c = 1. 

The vector and tensor perturbations have a much simpler form, which can be written in one concise equation as; 



with 2 



/ fin , -n , n(n-2) /iy/ T n 2 M\ Z v = 2,3 



W/r(r) = ^ I l« + »-!) + - ^^rr ^ = 1 ' 2,' .' .' .' ( L9 ) 



where \iy = 3 for the vector perturbations, while for the tensor perturbations \±t — — 1. Note that the mode I = 1 of 
the vector perturbation, which is absent from the above spectrum, represents a purely rotational mode of the BH. 

We shall mainly be interested in evaluating the scattering cross-section, which is related to the absorptions proba- 
bility, \A lP (u)\ 2 , by HI; 

(a \(n-l)/2 / J-. _i_ i \ 00 

= C [ -^ r ^-±-J E Dip\A ip {u)\ 2 =,5 + ^+^, (1.10) 



1 Note that for n = 2 the scalar (polar) perturbation agrees with the Zerilli equation |2l| . 

2 The vector (axial) perturbations for n = 2 reduces to the standard Regge- Wheeler equation |22 
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where P denotes each respective perturbation and the graviton normalization is |25| ; 



2 

C = ( n + 2 )(n-l) ' (LU) 



The degeneracy of each perturbation, Dip, is |2£j; 

(2l + n- l)(Z + n-2)! 



Dis 
Div 



(n-l)!Z! 
1(1 + n- 1)(2Z + n- 1)(7 + n - 3)! 
(n-2)!(i+l)! 



(n + l)(n - 2)(Z + n)(Z - 1)(2Z + n - 1)(Z + n - 3)! 
At = 2(n-l)!(J + l)! ' (L12) 

for a given angular momentum channel /, which is the spin-2 generalization of the result given in reference Q. Thus, 
only scalar and vector perturbations contribute for n = 2 and the partial sums effectively start from / = 2, even for 
the tensor perturbation. 



II. INTERMEDIATE ENERGY APPROACH 



To evaluate the absorption probability and hence the graviton emission rate in the intermediate energy regime we 
shall use the WKB approach of Iyer and Will |2jj . This was recently used in the higher-dimensional context by Berti 
et al. 13] (also see reference 0) to investigate the gravitational energy loss of high-energy particle collisions using 
a QNM analysis. In the following we shall work to lowest order in the generalised WKB method of reference |27j : 
however, to check the validity of the method we go up to second order to verify that the next order correction is small. 

As discussed in reference 13], in order to use the WKB method we must rewrite the perturbations in the (n + 2)- 
dimensional tortoise coordinate defined by; 



dr f(r) ' 

which implies that 



(2.1) 



'^Y.^LZA, (2-2) 



n — 1 ^— ' 



a- 



j=0 ~3 

where 

ay = (2M) 1 ^ n - 1 )e 27ri ^' /(n - 1) (j = 0, . . . , n - 2) . (2.3) 

The tortoise coordinate, r*, given above becomes quite complicated in more than four-dimensions and in our case it 
will be more convenient to work with the original coordinate r and use equation (|2.1|) to convert derivatives. Thus 
the perturbation equations reduce to the standard Schrodinger form; 

^ + Qp(r*)) $ = , Q P (r*) = u 2 - V P (n). (2.4) 

The subscript P denotes any one of the three possible perturbations. Note that the potential, equation (|1.3|) . is 
defined in terms of r and not r*, the tortoise coordinate. 

As discussed in reference 27], an adapted form of the WKB method can be employed to find the QNMs or the 
absorption probability, which we are primarily interested in, when the scattering takes place near the top of the 
potential barrier. In the following we shall use the same notation as reference 27]. The absorption probability, up to 
second order in the WKB expansion, is found to be; 

I^HI 2 = i + e J(^i/ 2 ) ' ( 2 - 5 ) 

where 

v + l/2 = i(2Q' Po )-^ 2 Q Po -A (2.6) 
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and 



A = i(2Q" )"V2 




Q 



Po 



32 I 3Q'/, o 



(v + l/2f 




30 I Q 



Po 



16 I 3Q% 



(2.7) 



these being parameters determined in reference 27], and where the subscript zero denotes the maximum of — Q(r*). 

By substituting equation (|2.7|) into equation (|2.6|) we can eliminate A to obtain an expression solely in terms of 
v + 1/2; however, because we wish to only know the size of the second order correction we can eliminate v + 1/2 to 
find A; 



A 



-L(l + 2AC±^l + A{A-B)c) , 



where 



and 



A = i(2Q%)- l ' 2 Q Po , 



(2, 



(2.9) 



-1/2 




c = «; )- 1 / 2 




(2.10) 

In general it is just a simple but tedious exercise in algebra to evaluate the second order correction. However, as we 
discuss in the next section, by considering the potential to leading order, O(l ), we can estimate the size of the second 
order correction without making any explicit calculations. 



III. ZEROTH-ORDER APPROXIMATION 



To illustrate the method as simply as possible let us work to zeroth order in the WKB expansion for the case of 
large /. However, in order to find a difference between the respective perturbations we should consider up to 0(1°) in 
the potential, which in this case is given by; 

V s(r) ^^{l 2 + in- 1)1 + + J(8 - 7n)(l - /)) , (3.1) 

for the scalar perturbation and 

Hvx(r) = ^ (> + (n - 1)/ + - ^n 2 (l - /)) , (3.2) 

for the vector and tensor potential, which is an exact expression to O(l ), see equation (|1.9|) . Importantly, we 
can immediately notice that in four-dimensions, n = 2, the scalar and vector potentials are identical to O(l ). The 
equivalence of the scalar and vector potentials for general / is well known in four-dimensions 28], however, as discussed 
in reference [20| |. this is not true for the case of general dimensions. 3 This will have important consequences for the 
graviton emission of a higher-dimensional BH, which we shall discuss later. For convenience let us write a general 
formula for all the perturbations, valid up to 0(1°); 

F P (r) = ^(/ 2 + ^-^(l-/)), where I 2 = I 2 + (n - 1)1 V/ = 2,3, ... . (3.3) 

In the previous equation we have used; 

P=^±, (3.4) 



3 To 0(1) all the perturbations become identical, regardless of the dimension. 
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and 

as = n(7n — 8), ay — 3n 2 , olt — —n 2 . (3.5) 

At this stage it will be convenient to change variables to z = c^r, for example see reference 29], and by doing so 
the WKB equation now becomes; 

jJ + Qp(*.)V = °> Qp(z.) = 1-Vp(z,), (3.6) 

where 

V, W _^(P + /J -^j£ T ), (3.7) 

and 

/CO = 1 - ; £ = . (3.8) 

The expression for the absorption probability, equation ((2.5)) . is defined explicitly in terms of derivatives of the 
potential with respect to the tortoise coordinate z* (or r*, where z* has been defined in an analogous way to r*); 
however, in terms of z, up to second order in derivatives of the potential, we find; 

dQ = dz dQ = ff z \ d Q = ^ ^\ d<3 , 3g x 

dz* dz* dz V z n_1 / dz 

% = >< ('<"§) = - 1^) <- + »> « f + o - ?M 0) • <-»> 

Note that at the maximum of the potential, zo, the first term in equation ((3.10)) is zero and thus simplifies our 
calculation. 

Hence, from equation ((2.6)) with A = 0, we are easily lead to the result; 

\^ip\ 2 = f T , (3.H) 

1 + e ¥ [f(zo)(P+(3-*tez i-")-z*] 

where we have written the general perturbation as; 

V2z 



,1+n 



K(l,z ,e) = e ^ y -2 n (1 + 2 n) a e 2 z 2 - 24 (f 2 + /?) z 2n + (2 + 3n + n 2 ) (4f 2 + a + 4/?) eV 

(3.12) 



and we have dropped the subscript P for clarity. Choosing the appropriate ap from equation ((3.5)) then determines 
the absorption probability for each respective perturbation mode. The location of the maximum, zo, for the potential 
— Qp(zo) can also be solved, see appendix A. This also depends on the type of perturbation. 

Before we proceed any further we should verify that the second order correction is small as compared to the zeroth 
order one. Given that Vp is of order 0(/ 2 ), see equation ((3.3)) . and that I is independent of r or r*, it is simple to see 
that; 

i(2Q' Po )~ 1/2 ~ 0{l~ l ) . (3.13) 

The terms in the square brackets for B and C, see equation ((2.10)) . are of O(l ) and hence (using the relation above) 
B and C contribute to order 0(/ _1 ). Furthermore the term A is of order 0(1). Thus the second order contribution, A, 
see equation ((2.8)) . is of order (9(/ _1 ) and is therefore negligible in the limit I ^> 1. This is also true for the high-energy 
limit e ^> 1, as well as for intermediate energies e ~ l. 4 



4 The WKB method breaks down for the low-energy case e <C 1. 
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FIG. 1: Black hole absorption probability, |*4zp| 2 , as a function of angular momentum / (extended to real numbers) for e — 1. 
Blue (solid), red (dashed), black (dot-dashed) and purple (dotted) curves correspond to n — 2,3,4 and 5 respectively. Note 
that the scalar and vector perturbations are only valid from / = 2 onwards, however, for comparison with the tensor case we 
start at / = 1. 



IV. RESULTS 



In this section we shall present our results as based on the analytic expressions derived in the previous section. 
Furthermore we shall also demonstrate that the geometric optics limit is reproduced for large values of e. Note that 
for completeness we shall also derive this expression. 

We have presented our first set of results in Fig. [I] where we have plotted the absorption probability as a function of 
I for each of the perturbations. We have also considered scenarios involving different dimensions, n, for e ~ 1. Firstly, 
the scalar mode has the largest contribution, followed by the vector and then the tensor perturbations; however, 
this does depend on the value of e. Secondly, we see that the absorption probability is larger for smaller n. In fact 
numerical plots for different e show that |v4jp| 2 saturates to unity for larger and larger I as e increases. This implies 
that larger I is required in the angular momentum sum for the cross-section, when obtaining the geometric optics 
limit, see below. 

Turning to the cross-section, as given by 5 ; 



, M = c ( 47r) ("-^ r (H±l) (™) - 1 V V . ^ r , (4.1) 

where Dip is given in equation (|1.12|) and remembering that a = a p. It should be noted that though we are primarily 
interested in the intermediate energy regime, e ~ 1, with large angular momentum, I ^> 1, the above expression can 
also be applied to the case where I > s 1 /^ -1 ) ~ 1. Note also that although the next order correction in the WKB 
expansion becomes larger for I ~ 0(1) we can still extrapolate the results hoping the errors are not too large. In 



5 The cross-section effectively starts from / = 2 for the tensor perturbations, like for the other perturbations, because of the form of the 
degeneracy Z^t- 
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FIG. 2: Black hole cross-sections as a function of e, in units of (2M) , for each gravitational perturbation and the total 
one. The horizontal lines correspond to the classical cross- section. Blue (solid), red (dashed) and black (dot-dashed) curves 
correspond to n — 2, 3 and 4 respectively. Note that there is no tensor perturbation for the n — 2 case. 



Fig. [2 we plot the cross-section, equation (|4.1[h as a function of e and to ensure convergence in the partial wave sums 
we sum up to / 

max ^ SSmax' The total cross-section tends to the classical one for large £ as we shall now explain. 
In the intermediate energy approach (for large /) the absorption probabilities satisfy Ais = A\y — Ait, see Fig. [I] 
and the degeneracies satisfy [25l | : 

Dis ^ C {D ls + D lv +D iT ) . (4.2) 

Furthermore, the classical cross-section corresponds to the high-energy limit, ^/O- 1 ) ^> / ^> 1, which then implies 
the absorption coefficients Aip ~ 0(1) and in this limit the sum over / has a cut-off at / « boo (for example, see 
reference |23) where b is the critical radius (or impact parameter) at which the BH ceases to absorb radiation. This 
is given by equation (12) of reference 31], for a massless/relativistic particle as; 



b ={—) {—i) (2M) ^- (43) 

Given the cut-off in the mode sum we have; 

I 

for large a;, and substituting this into equation (|4.1[) . using (|4.2|) . leads to the classical cross-section; 

cj n V 2 / ™ ! nT{n/2) \ 2 ) \n - 1 J K J ' v; 

where in the second step we used the gamma duplication formula. This agrees with the standard result in (2 + n)~ 
dimensions, compare with reference 3 (after making the substitution r^f 1 = 2M). 
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As is well known, the high-energy cross-section is independent of the particle species, and likewise we see that it is 
independent of the graviton cross-section. These are represented by horizontal lines in the bottom right plot in Fig. [2| 
and we see that our analytic results correctly reproduce the high-energy limit e ^> 1 (note that for higher-dimensions 
larger values of e are required to obtain the geometric optics limit 6 ). 

V. CONCLUSION 

In conclusion we have shown in this paper that the WKB method of Iyer and Will |23] can be applied to the case of 
graviton emission from a Schwarzschild BH. Indeed our results reproduce the classical cross-section in the high-energy 
limit. We have also presented new results for the intermediate energy regime. 

In the low-energy limit, Moo <C 1, our method breaks down as can be seen from Fig. 2, where the cross-section starts 
to diverge. It is well known that the cross-section is proportional to oj 2l + n in (2 + n) -dimensions Q, which essentially 
corresponds to s- wave scattering for a spin-zero field, as the lowest / modes dominate the cross-section. Indeed it 
is straightforward to obtain the low-energy cross-section for the vector and tensor perturbations by employing the 
standard tec hniq ue of matching the near horizon and far field solutions, for example using the techniques discussed 
in references |32l l33l l34|. However, there is a subtlety with the higher-dimensional Zerilli (scalar) equation: In four- 
dimensions the Regge- Wheeler (vector) solution is usually used to find the absorption probability, which is equivalent 
to that for the Zerilli equation, as the two solutions are identical in four-dimensions 28], but as pointed out in reference 
[20j | . in higher dimensions no such relation exists. Thus, if we attempt to use the standard technique of matching 
the near horizon and far field solutions directly for the Zerilli equation it appears that numerical methods seem more 
amenable. Regardless, the focus of this current work has been the intermediate energy regime given that this is where 
we expect graviton emission to be most interesting for a rapidly rotating BH. 

Although in this article we focused on the static Schwarzschild BH, we can also apply our method to the case of 
spin-zero field emission from a Kerr BH 35] . Note that though the solution to the graviton perturbations for a higher- 
dimensional Kerr BH have not yet been found, we are encouraged by the fact that there are some similarities between 
the total graviton and spin-0 cross-sections, see Appendix B. As can be seen from Fig. [31 the total cross-section, which 
is the sum of the perturbations (scalar, vector and tensor) is of the same order of magnitude as the spin-0 case, i.e. 
0(<js + &v + &t) ~ 0(gb), and such an approximation may also be valid for the rotating case. 

The rotating case is particularly interesting due to the phenomenon of super-radiance (for example see reference 
[34)), where the absorption probability becomes negative. Super-radiance has been discussed in the higher-dimensional 
context recently in references 0, 0, 0, 0. However, the highly rotating case has only been considered 
using numerical techniques; whereas we expect that our approach allows for an analytic expression for such a case. 
Furthermore, our approach should also allow one to evaluate the QNMs analytically for the rapidly rotating case, as 
has recently been done in the limit a — > 00 |36| . In the four-dimensional case the rotation parameter a = 1 is bounded 
to a max = 1, for example see reference [28|. This fact has meant that there have been no prior investigations of the 
QNMs of a Kerr BH with a large rotation parameter. However, in (2 + n)-dimensions there are (n + l)/2 rotation 
parameters, where for n > 3 the rotation parameters are unbounded [24(. 7 Indeed, as we mentioned, even after the 
spin-down phase in some models the BH remains rotating 

Another interesting case for which the WKB method can be applied is that for emission from a charged (Reisner- 
Nordstrom) BH. Recently the graviton perturbations have been found in reference I37J (along with their QNMs [Tl^. 
However, in this case the solutions are complicated by the fact that there are now solutions for the electromagnetic 
field itself. We hope to present results for this model in the near future. 
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The geometric optics limit in equation 14.51 has a maximum at n « 7, which implies that for n > 7 (i.e. for greater than 9-dimensions) 
the classical cross-section starts to decrease. 

If we assume that the BH is created on a thin brane then there is only one rotation parameter, for example see reference 
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Note added: After the completion of this article, similar work appeared in 25], which correctly pointed out that 
the degeneracy factor for the graviton cross-section should be that for a spin-2 field [2(|, which is different to the 
spin-0 case, for n > 2. Our main results do not change; however, only the total cross-section now coincides with the 
classical one, due to the normalisation. This actually improves our attempt to model the spin-2 perturbations by a 
spin-0 field, as is discussed in Appendix B. 



APPENDIX A: LOCATION OF THE MAXIMUM, Q Pq 

In order to find the maximum of the potential it is convenient to work with the coordinate; 

* = W = ^ (AD 
rather than r* or r. That is we wish to find the roots of; 

ar* ar* dr dx 

Here we should mention that in our case, even though all the expressions depend on Qp, the potential can be written 
in the form Qp = uo 2 — Vp and thus, for a given energy, cj, the maximum of — Qp corresponds to the maximum of Vp. 
Therefore, apart from the solutions at the horizon (x = 1) and at infinity (x = 0), we merely require the roots of; 

Working to 0(1) in the potential allows the maximum to be easily found; 

2 



n+ 1 ' 



(A4) 



which is a result valid for all three perturbations. However, as we have previously discussed, in order to see a difference 
between the perturbations we must work to order O(l ). Hence, for the gravitational perturbations we obtain; 

(1 + n) (4/ 2 + a + 4/^ - \ -Una ( P + (3) + (1 + nf (ll 2 + a + 4(3) 2 

— - 1 — — -• 

where, as before, we are using the shorthand notation I 2 = I 2 + (n — 1)1 and the coefficients a and (3 (ignoring the 
subscripts) as defined in equation (|3.5j) . Note that if we expand the square root in powers of large I we obtain equation 

(1X1 . 

The location of the maximum in terms of the coordinates r and z is then simply given by; 
which is independent of the perturbation up to 0(1), see equation (|A4|) . 



APPENDIX B: SPIN-ZERO EMISSION 



In this appendix we briefly discuss the case of emission of a bulk scalar (spin-zero 
potential is equivalent to that for the tensor perturbations. Following reference 
the scalar radial equation, g^ v V^Vz, <j)(x) = 0, is found to be; 



oH 
32] 



field. As we shall see the effective 
after a separation of variables, 



f(r) d 



f(r) r 



n dR(r) 
dr 



fir) 



1(1 



1) 



R(r) = . 



(Bl) 
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FIG. 3: Black hole cross-section for the graviton compared to a spin-zero field as a function of e, in units of (2M) n ~ 1 . The 
horizontal lines correspond to the classical cross-section. Blue (solid), red (dashed) and black (dot-dashed) curves correspond 
to n = 2, 3 and 4 respectively. 



This equation can be written in the WKB form by applying the tortoise coordinate defined in equation (|2.1|) 8 and 
rescaling the radial solution by R(r*) = r _n / 2( £(r*). A short calculation leads to; 



Z(Z + n-l) + ^(n-2) + ^(l-/) 



*(r.) = 0, 



Ib = 0,1,2,. 



(B2) 



One can then easily verify that this equation is identical to that for the tensor perturbations of the graviton, see 
equation (|1.9|) : however, the degeneracy factor, Dip, will be different. The only difference between the two being that 
the sum over / starts from / = for a spin-zero field, not 1 — 2. 

In Fig. [31 we compare the spin-zero cross-section, ag, with the total graviton cross-section. Although the classical 
cross-sections agree in the high-energy limit (asymptotically), at intermediate energies they do not agree due to 
differences in the partial wave sums. However, at lower energies we see that by taking the spin-0 cross-section from 
/ = 1 we obtain the best agreement with that for the graviton, see Fig. 03 Note that for further comparison we also 
plot the spin-0 cross-section with the partial sum starting from 1 = 2. 



In reference lo'J a slightly different choice of tortoise coordinate is made. 
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